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Abstract 

A  dynamic  model  for  a  stand-alone  renewable  energy  system  with  hydrogen  storage  (RESHS)  is  developed.  In  this  system,  surplus  energy 
available  from  a  photovoltaic  array  and  a  wind  turbine  generator  is  stored  in  the  form  of  hydrogen,  produced  via  an  electrolyzer.  When  the 
energy  production  from  the  wind  turbine  and  the  photovoltaic  array  is  not  enough  to  meet  the  load  demand,  the  stored  hydrogen  can  then 
be  converted  by  a  fuel  cell  to  produce  electricity.  In  this  system,  batteries  are  used  as  energy  buffers  or  for  short  time  storage.  To  study  the 
behavior  of  such  a  system,  a  complete  model  is  developed  by  integrating  individual  sub-models  of  the  fuel  cell,  the  electrolyzer,  the  power 
conditioning  units,  the  hydrogen  storage  system,  and  the  batteries  (used  as  an  energy  buffer).  The  sub-models  are  valid  for  transient  and  steady 
state  analysis  as  a  function  of  voltage,  current,  and  temperature.  A  comparison  between  experimental  measurements  and  simulation  results  is 
given.  The  model  is  useful  for  building  effective  algorithms  for  the  management,  control  and  optimization  of  stand-alone  RESHSs. 

©  2004  Elsevier  B  .V.  All  rights  reserved. 
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1.  Introduction 

The  use  of  a  stand-alone  renewable  energy  system  (i.e.  one 
using  wind  and  photovoltaic  energy)  in  remote  areas  requires 
an  energy  storage  device  to  smooth  out  the  intermittent  power 
input  from  these  sources.  Recent  system  designs  [1-4]  rely 
on  batteries  for  short-term  energy  storage,  while  hydrogen  is 
used  for  long-term  energy  storage.  In  these  systems,  the  hy¬ 
drogen  (EL)  has  been  produced  through  an  electrolyzer  pow¬ 
ered  by  the  surplus  energy  available  from  the  primary  sources 
(wind  turbine  and  photovoltaic  array).  When  the  input  power 
is  insufficient  to  feed  the  RESHS  load,  previously  stored  hy¬ 
drogen  is  reconverted  through  a  fuel  cell  (FC)  to  produce  the 
required  electricity.  The  design,  management  and  optimiza¬ 
tion  of  such  a  system  require  a  useful  model. 
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We  present  a  model  to  describe  the  dynamics  of  an 
RESHS.  It  integrates  sub-models  of  the  electrolyzer,  the  fuel 
cell,  the  batteries,  the  power  interfaces  (buck  and  boost  con¬ 
verters)  and  the  storage  system.  Interdependency  issues  (hy¬ 
drogen  consumption  cannot  exceed  production)  are  taken  into 
account.  Special  attention  is  given  to  the  characterization  of 
the  system’s  major  components  in  the  transient  state,  and  we 
use  simple  and  realistic  assumptions  to  describe  the  behav¬ 
ior  for  short-  and  long-term  operation  of  the  RESHS.  Most 
of  the  sub-models  are  specified  by  the  component’s  polar¬ 
ization  curves  characteristics  (current- voltage-temperature). 
The  model  is  validated  by  comparing  its  output  to  that  of 
the  Hydrogen  Research  Institute’s  (HRI)  renewable  energy 
system  test  bench,  which  is  completely  described  in  [4,6] 
and  whose  configuration  and  specifications  are  given,  respec¬ 
tively,  in  Fig.  1  and  Table  1.  A  scenario  built  with  realistic 
residential  power  consumption  needs  and  typical  power  pro¬ 
duction  by  wind  turbine  (WT),  and  photovoltaic  (PV)  array 
is  also  simulated  and  analyzed. 
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Fig.  1.  Block  diagram  of  the  renewable  energy  system  test  bench  of  the  HRI. 


2.  Modeling  of  the  components 

Generally,  a  RESHS  is  designed  for  a  nominal  dc  bus  volt¬ 
age,  which,  in  the  case  of  the  HRI  test  bench,  is  about  48  V. 
However,  the  real  voltage  on  the  dc  bus  depends  on  the  oper¬ 
ating  conditions  of  the  system.  When  the  energy  production 
exceeds  what  is  needed  and  the  battery  (short-term  energy 
storage  device)  is  being  charged,  the  input  power  devices 
tend  to  impose  their  output  voltage  on  the  dc  bus.  Wind  gusts 
can,  for  example,  increase  that  bus  voltage  from  48  V  to  56  V 
in  a  fraction  of  a  second.  Similarly,  when  input  energy  pro¬ 
duction  is  below  what  is  needed  and  the  load  draws  on  the 
battery,  it  is  the  battery  that  will  impose  its  voltage  on  the 
dc  bus.  This  variability  of  the  bus  voltage  is  a  major  control 
problem,  as  quite  clearly  this  voltage  cannot  be  considered 
as  a  reliable  variable  to  describe  the  evolution  of  the  state  of 
the  RESHS.  Instead,  it  is  the  battery  energy  that  will  be  used 
as  a  system-controlling  variable  (see  Section  2.1). 


Table  1 


RE  test  bench  technical  specifications 


Components 

Type 

Power  (kW) 

Voltage  (V) 

Photovoltaic  (PV) 

1 

48 

Wind  turbine  gener- 

10 

48 

ator  and  regulator 

Electrolyzer 

Alkaline 

5 

26-48 

Buck  converter 

Multiphase  PWM 

5 

26-48 

Fuel  cell  (FC) 

PEM 

5 

24 

Boost  converter 

Multiphase  PWM 

5 

24-48 

Inverter 

5 

110  ac 

Load 

0-5 

110  ac 

Capacity  (kWh) 

Batteries 

Lead-acid 

10.5-55 

48 

Storage  H2 

125 

PWM:  pulse  width  modulation;  PEM:  proton  exchange  membrane. 


In  this  paper,  most  of  the  models  are  described  as  functions 
of  time,  current,  voltage,  and  temperature.  For  simulation 
purposes,  the  input  signals  are  the  wind  generator  rectifier 
output  current  (Avt),  the  PV  array  regulator  output  current 
(/py),  and  the  load  current  (/l).  Due  to  the  intermittent  nature 
of  the  renewable  energy  sources,  sampled  signals  will  be  used 
to  represent  all  of  them.  This  way,  any  energy  production  and 
load  profile  can  be  modeled  at  will.  In  the  following  sections, 
the  models  of  the  sub-units  are  presented  in  the  order  in  which 
they  are  traversed  by  the  energy  flux:  battery,  buck  converter, 
electrolyzer,  boost  converter,  fuel  cell,  and  hydrogen  storage. 


2.7.  Battery  model 


The  battery  is  the  main  component  on  the  dc  bus,  and  plays 
the  role  of  an  energy  buffer  to  handle  current  spikes  and  for 
short-term  energy  storage.  Different  models  for  batteries  are 
available,  in  particular  those  suitable  for  electrical  vehicle 
applications  [5,11,12,15].  For  stationary  applications,  such 
as  the  RESHS,  the  models  described  in  [2]  use  many  experi¬ 
mental  parameters  that  cannot  be  estimated  easily,  such  as  the 
overcharge  effect  (though  in  a  properly-controlled  RESHS, 
this  effect  does  not  happen,  and  hence  is  not  included  in 
the  model).  The  main  parameters,  which  determine  the  bat¬ 
tery’s  performance,  are  its  internal  resistance,  the  polariza¬ 
tion  effect,  and  the  long-term  self-discharge  rate.  This  self¬ 
discharge  rate  is  difficult  to  estimate,  and  is  itself  subject  to 
a  number  of  factors,  such  as  the  operating  temperature,  the 
number  of  operation  cycles,  and  the  materials  and  technology 
used  in  its  manufacture  [9,14]. 

The  battery  voltage  7/b  it),  which  takes  these  three  param¬ 
eters  into  account  is  given  by, 


UB(t)  =  (1  +  at)UB,  o  +  Ri(t)I(t)  +  K{QR(t)  (1) 

where  a  is  the  self-discharge  rate  (s-1);  7/b,0  is  the  open 
circuit  voltage  (V)  at  t  =  0;  Rft)  is  the  internal  resistance  (£2), 
K[  is  the  polarization  coefficient  (£2h_1),  and  (9r(0  is  the 
rate  of  accumulated  ampere  hours.  If  I(t)  >  0  then  the  battery 
is  charging;  if  I(t)  <  0  then  the  battery  is  discharging.  The 
battery  energy  is  then, 

W(t)  =  Wo  +  f  Pin(t')dt'  (2) 

J  0 

where  P-mif)  =  UB(t)I(t)  is  the  input  power  to  the  battery  and 
Wo  is  the  battery’s  initial  energy.  As  we  will  see  later,  the 
decision  algorithm  (as  to  whether  electrolyzer  or  fuel  cell  are 
to  be  activated  to  rebalance  the  battery  energy)  will  depend 
on  the  battery’s  state  of  charge  (SOC),  defined  by, 


SOC(0  = 


W{t) 

Wmax 


where  Wmax  is  the  maximum  battery  energy  without  over¬ 
charge. 
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Fig.  2.  Diagram  of  the  identification  process  of  the  buck  converter  parame¬ 
ters. 

2.2.  Buck  converter 


When  there  is  an  excess  of  electrical  energy  in  the  system, 
that  excess  is  channeled  to  the  electrolyzer  to  produce  hydro¬ 
gen.  To  control  this  hydrogen  production,  a  buck  converter 
designed  at  HRI  controls  the  input  current  to  the  electrolyzer 
cells.  This  buck  converter  is  a  dc  voltage  reducer  designed 
to  maximize  the  power  transfer  from  the  dc  bus  to  the  elec¬ 
trolyzer  cells.  It  uses  the  multiphase  technique  to  generate 
pulse-width  modulated  (PWM)  signals  [6].  Its  models  are 
expressed  by  Eqs.  (4)  and  (5),  and  they  give  the  voltage  and 
current  applied  to  the  electrolyzer  cells.  Considering  the  com¬ 
plexity  of  describing  these  circuits,  an  autoregressive  moving 
average  (ARM A)  model  (Fig.  2)  [16]  for  parameter  identifica¬ 
tion  is  used  to  minimize  the  average  quadratic  error  between 
the  real  output  of  the  buck  converter  and  the  model’s.  During 
this  process,  the  time  scale  is  set  to  1  s  and  the  buck  converter 
operating  temperature  to  25  °C. 

The  following  equation  (based  on  the  ARMA  model)  gives 
the  relation  between  the  buck  output  voltage  and  the  dc  bus 
voltage  (transfer  function), 


^Bu,Out(ft)  =  Ub  in) 


(^Bu,0  T  Bbu,1Z  1)  Dgu ifl) 
Abu,0  +  Abu,1  Z~l 


Table  2 

Parameter  values 


Component 

Parameters 

Values 

Battery 

Ri  («) 

0.076 

Ki 

~0 

Qfc  (Ah) 

880 

Qo  (Ah) 

880 

£o(V) 

48 

Boost  converter 

Abo,o 

1 

Abo,i 

-0.2045 

^Bo,o 

1.236 

^Bo,i 

-0.4275 

Buck  converter 

Abu,o 

1 

Abu,i 

-1.485 

*71, el 

0.7 

^Bu,o 

0.057 

^Bu,i 

-0.082 

Electrolyzer 

1,0 

22.25 

Cl 

-0.1765 

C2 

5.5015 

*71, FC 

0.45 

4i,o 

0.1341 

^el 

-3.3189 

A Cell, el 

24 

PEMFC 

Cfc 

33.18 

Ei 

-0.013 

e2 

-1.57 

4c, o 

8.798 

^FC 

-2.04 

A Cell,FC 

35 

verter  is 


^Bu,In(/0 


U Bu,  Out  (ft )  f Bu,  Out  (ft ) 
^7Buf^Bu,In(^) 


where  /Bu,Out  is  the  input  current  to  the  electrolyzer  cells, 
and  is  determined  in  the  next  section.  The  coefficients  may 
vary  slightly  from  one  operating  point  to  another,  but  on  the 
whole,  the  model  gives  a  good  account  of  the  dynamics  of 
the  system. 


2.3.  Electrolyzer 


As  the  electrolyzer  time  response  is  slow  [10]  compared  to 
the  modeling  sampling  time  (1  s),  the  output  voltage  is  given 

by, 

//ei(0\  Rcl 

UQl(t)  =  UQ  i,o  +  CM)  +  C2ln  +  ^4l(0 

\  4l,0  J  f  elw 

(6) 


where  Abu,o,  Abu,i,  #Bu,o,  and  Z?bu,i  are  parameters  which 
have  to  be  determined.  Their  values,  for  our  case,  are  given 
in  Table  2.  Ub  is  the  dc  bus  voltage;  Dbu  is  the  duty  cycle 
and  f/Bu,Out  is  the  buck  converter  output  voltage  (and  applied 
to  the  electrolyzer  cells).  The  buck  converter  input  voltage 
(^Bujn)  is  equal  to  the  dc  bus  voltage  because  of  its  direct 
connection  to  the  dc  bus.  Taking  into  account  the  buck  power 
efficiency  (r/Bu),  the  input  current  (/Bujn)  to  the  buck  con- 


where  t/d,o  (V),  C,  (V  C'-'),  C2  (V  °C-1),  /ei,o  (A)  and  Rel 
(£2°C-  )  are  parameters  to  be  determined  experimentally. 
For  the  HRI  test  bench,  their  values  are  given  in  Table  2  for 
an  operating  temperature  Te\(t)  between  22  °C  and  52  °C.  The 
first  two  terms  of  Eq.  (6)  represent  the  theoretical  potential 
of  an  ideal  cell.  The  third  term  gives  the  activation  potential, 
while  the  last  one  represents  resistance  (ohmic)  losses.  These 
parameter  values  vary  from  one  electrolyzer  to  another.  Con- 
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sidering  that  the  output  of  the  buck  converter  is  connected 
directly  to  the  input  of  the  electrolyzer  cells,  t/Bu,Out(0  and 
^Bu,Out(0  are  equivalent  to  UQ\(t)  and  IQ\(t),  respectively.  The 
hydrogen  production  rate  Vei(0  is  given  by, 


Vel  =  Ncell,el 


?7l,el4l(0 

CH, 


where  771^1  is  the  electrolyzer  utilization  factor;  A^cellel  repre¬ 
sents  the  number  of  cells,  and  Ch2  is  a  conversion  coefficient 
(of  value  2.39  A  hi-1  of  hydrogen).  t]\,q\  depends  on  the  cell 
temperature.  For  alkaline  electrolyzers,  some  studies  [7,8,13] 
show  that  771^1  is  between  0.6  and  0.75.  A  more  rigorous  es¬ 
timate  can  be  obtained  by  direct  measurement  of  the  produc¬ 
tion.  The  hydrogen  produced  by  the  electrolyzer  is  stored,  and 
can  be  used  later  by  the  FC.  The  storage  sub- system  model 
is  given  in  Section  2.5. 


measured  on  the  IRH  test  bench  are  given  in  Table  2  for  tem¬ 
peratures  ranging  from  24  °C  to  72  °C.  7pc(0  is  the  operating 
temperature  of  the  cells,  which  changes  during  the  cell’s  op¬ 
eration.  The  first  two  terms  of  Eq.  (10)  represent  the  open 
loop  potential,  the  third  term  corresponds  to  the  activation 
potential,  and  the  last  represents  resistance  (ohmic)  losses. 
These  parameters  vary  widely  from  one  fuel  cell  to  another. 
The  FC’s  hydrogen  consumption  rate  Ffc(0  is  given  by, 


VFC  =  Acell,FC 


?7l,FC^FC(0 

Cu2 


(11) 


where  rji FC  is  the  utilization  factor  of  the  FC,  and  A^celLFC 
represents  the  number  of  cells.  The  literature  [7,8]  shows 
that  rji fc  is  between  0.3  and  0.6  for  PEM  fuel  cells.  A  more 
accurate  estimate  can  be  obtained  by  direct  measurement. 


2.4.  Boost  converter 


2.6.  Hydrogen  storage  sub-system 


The  nominal  output  voltage  of  HRFs  fuel  cell  is  around 
24  V,  whereas  the  dc  bus  is  designed  for  48  V.  Hence  a  power 
interface  is  necessary  between  the  two.  In  order  to  maximize 
the  power  transfer  between  the  FC  and  the  dc  bus  the  boost 
converter  is  designed  with  the  same  multiphase  switching 
technique  and  pulse- width  modulated  signals,  as  was  used 
for  the  buck  converter.  Using  the  same  method  as  described 
in  Section  2.2,  the  relation  between  the  duty  cycle  D^0  and 
the  input  current  /Bojn(0  of  the  boost  converter  is  given  by, 


7Bo,In(ft)  =  7pc,Max 


(#Bo,0  +  #Bo,U  l)  DBo(n) 
Abo,0  +  Abo  ,lZ~l 


where /pc.Max  is  the  maximum  output  current  oftheFC;  Ab0,o> 
Abo,i,  #Bo,o>  and  #bo,i  are  parameters  to  be  determined  (see 
Table  2).  The  output  current  (/Bo,Out)  of  the  boost  converter 
is  obtained  from  the  boost  power  efficiency  (ijbo) 


^Bo,Out(^)  —  7?Bo 


Thin) 


where  I/fc  is  the  FC  output  voltage,  and  ijBo  is  determined 
by  direct  measurement  (and  found  to  be  >95%). 


The  hydrogen  produced  by  the  electrolyzer  is  stored  in  a 
pressurized  tank.  Although  the  actual  renewable  energy  sys¬ 
tem  at  HRI  has  a  compressor  and  the  hydrogen  can  be  stored 
in  a  tank  at  up  to  150psi,  the  equations  used  will  assume 
hydrogen  storage  at  normal  pressures  and  temperatures.  The 
power  balance  for  the  storage  function  is  given  by 


dWk2(0 

d  t 


=  Pclit) 


Pfc  (0 


(12) 


where  Wu2  ( t )  is  the  energy  present  in  the  tank  in  the  form  of 
hydrogen,  Pfc(0  and  Pe\(t)  are,  respectively,  the  FC  power 
consumption  (related  to  the  H2  consumption  rate)  and  the 
electrolyzer  power  production  (related  to  the  H2  production 
rate).  The  solution  of  Eq.  (12)  in  the  Laplace  domain  is  given 
by 


Wn2(s) 


Pq \(s)  ~  Pfc(s) 

- ; - +  wh2,0 


(13) 


with 

n  ,  s  %\(s)  VfcCO 

Pd(s)  =  — — — AH;  PFC(s)  =  ——AH 

Vj  Vj 


2.5.  Proton  exchange  membrane  fuel  cell 


According  to  references  [6-8],  the  proton  exchange  mem¬ 
brane  (PEM)  FC  reaction  time  constants  during  transient  pe¬ 
riods  (<50  ms)  are  smaller  than  the  simulation  step  time  (Is). 
Hence,  the  dynamic  behavior  of  the  FC  is  given  by  its  polar¬ 
ization  curve 


Ufc (t)  =  Ufc,o  +  E\  7fc(0  +  E2  In 


+ 


^FC 

Tpc(t) 


Ifc  (0 


(10) 


where  UFC, o  (V),  E i  (V0^1),  E2  (V0^1),  /FC,o  (A),  and 
Rfc  (£2  °C_1)  are  experimental  parameters  whose  values  as 


where  Vei  and  Vfc  are,  respectively,  the  production  and  the 
consumption  rates  of  H2.  Wu2  0  is  the  initial  hydrogen  en¬ 
ergy  stored.  Under  normal  temperature  and  pressure  condi¬ 
tions  [9],  Vt  is  equal  to  22.41  mol-1.  AH  is  the  enthalpy  of 
hydrogen  (“high  heating  value”)  (AH  =  286  kJ  mol-1). 


3.  RESHS  simulation 

A  model  of  a  RESHS  has  been  developed  based  on  the 
above  equations  for  its  sub-units.  The  simulation  block  dia¬ 
gram  is  given  in  Fig.  3.  As  mentioned  earlier,  the  simulation 
and  sampling  time  step  is  taken  as  1  s,  and  A  is  the  total  num¬ 
ber  of  steps  for  one  complete  operation  period.  The  principal 
results  of  interest  are  the  output  currents  of  the  buck  converter 
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Fig.  3.  RESHS  simulation  block  diagram. 

and  of  the  boost  converter,  and  the  energy  stored  in  the  bat¬ 
tery.  The  boost  converter  output  current  depends  of  course  on 
the  details  of  the  models  of  the  boost  converter  and  of  the  FC. 
Similarly,  the  electrolyzer  current  is  closely  associated  with 
the  details  of  the  buck  converter  and  of  the  electrolyzer.  In 
this  diagram,  Q(n)  represents  the  accumulated  ampere  hours 
in  batteries. 


4.  Experimental  and  performance  analysis 

The  models  of  the  sub-units  described  previously  were 
used  to  put  the  RESHS  system  through  a  typical  run  (the 
values  of  the  parameters  of  the  various  sub-units  are  found 
in  Table  2).  In  order  to  reflect  the  variable  nature  of  the  ac¬ 
tual  WT  and  PV  array  power  variations,  the  simulation  inputs 
were  subjected  to  relatively  large  ripples,  and  the  simulation 
outputs  were  monitored  closely  to  determine  the  stability  of 
the  model.  The  simulation  inputs  are  the  WT  generator  cur¬ 


rent  at  the  dc  bus  (/wt)>  the  PV  array  regulator  current  (/py), 
the  load  current  (7l),  the  electrolyzer  cell  temperature  (rei) 
(which  is  relatively  constant:  25  °C)  and  the  FC  cell  temper¬ 
ature  (7pc),  which  varies  linearly  from  25  °C  to  35  °C  with 
a  slope  of  3.3  °Cs-1.  The  simulation  outputs  are  (as  men¬ 
tioned  previously)  the  buck  converter  output  current  (same 
as  the  electrolyzer  input  current)  (7ei),  the  boost  converter 
output  current  (7B0,Out)>  and  the  battery  energy  (W(t)).  All 
these  signals  were  also  measured  on  the  actual  HRI  operat¬ 
ing  test  bench  to  permit  comparison  of  the  simulation  output 
with  reality. 

The  electrolyzer  and  FC  regulators  work  so  as  to  bring  the 
battery  energy  W(t)  back  towards  its  reference  value,  Wref, 
whenever  Wit)  gets  beyond  a  pre-determined  range  bounded 
by  Whi  and  Wiow  bracketing  Wref.  The  electrolyzer,  which 
lowers  W(t)  by  converting  the  excess  battery  energy  into  hy¬ 
drogen  (i.e.  when  the  input  energy  is  greater  than  what  is 
needed  by  the  load),  is  powered  on  when  W(t)  >  Whi,  and  is 
stopped  when  Wref  is  reached.  The  FC,  which  raises  W(t)  by 
converting  hydrogen  back  into  battery  energy  (i.e.  when  the 
input  energy  is  smaller  than  what  is  needed  by  the  load),  is 
powered  on  when  W(t)  <  W\ow ,  and  is  stopped  again  when 
Wref  is  reached.  The  system  starts  with  an  initial  energy  Wo 
>  Whi  (Wo  is  estimated  to  42,240 Wh),  usually  at  Whi.  In 
the  first  simulation  Wref  =  42,218  Wh,  and  Whi  and  Wiow 
are  at  42,240  Wh  (=Wref  +  22  Wh)  and  42,1 96  Wh  (=Wref 
—  22  Wh).  The  control  range  thus  extends  to  ±22  Wh  around 
Wref. 

4.1.  Input  signals 

Figs.  4  and  5  present  typical  WT  output  and  load  current 
(two  of  the  simulation  input  variables).  It  can  be  noticed  that 
the  system  operates  during  7  min.  There  was  no  significant 
power  available  from  PV  array  when  running  this  experience. 


4.2.  Electrolyzer 

The  electrolyzer,  which  is  fed  from  the  buck  converter,  is 
set  to  operate  at  a  nominal  input  power  of  1500  W  when  it  is 


Fig.  4.  Current  /wt(0  from  the  wind  turbine  through  the  rectifier  module. 
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Fig.  5.  Profile  of  the  load  current. 

running.  When  the  simulation  starts,  the  electrolyzer  is  oper¬ 
ating  according  to  the  control  strategy  described  above.  Fig.  6 
compares  the  experimental  and  the  simulation  results.  The 
difference  between  simulation  and  measurements  is  around 
8.5%.  This  result  shows  that  the  models  of  the  buck  converter 
and  the  electrolyzer  are  relatively  realistic.  The  steady- state 
difference  between  the  two  (at  t  >  65  s)  is  due  to  the  inaccu¬ 
rate  estimation  of  the  initial  energy  of  the  battery  (Wo)  in  Eq. 
(2)  (it  turns  out  that  the  model  of  the  energy  buffer  is  itself 
dependent  on  the  initial  energy  stored  in  the  battery).  In  spite 
of  this,  the  start-stop  sequences  of  the  electrolyzer  are  the 
same  for  the  simulation  and  the  experimental  measurements. 

4.3.  Fuel  cell 

For  the  FC  and  the  boost  converter,  the  simulated  and 
experimental  (measured)  currents  are  perfectly  matched 
(Fig.  7).  The  deviation  between  the  two  is  less  than  2%. 
It  turns  out  that  the  boost  converter  output  current  is 
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Fig.  7.  Comparison  of  the  simulated  and  measured  current  of  the  boost 
converter  output. 


not  very  sensitive  to  reasonable  variations  of  the  dc  bus 
voltage. 


4.4.  Buffer  energy 

The  overall  quality  of  the  simulation  of  the  operation  of 
the  RESHS  is  judged  by  examining  the  time-evolution  of  the 
buffer  energy  ( W(t )).  Fig.  8  gives  the  measured  and  simulated 
results.  From  f=0stof=65s,  W(t)  decreases  quickly,  because 
(1)  insufficient  energy  is  available  from  the  wind  turbine  and 
the  PV  array  (see  Figs.  4  and  5),  and  (2)  the  electrolyzer  is 
running  and  withdrawing  energy  from  the  battery  (Fig.  6). 

When  W(t )  reaches  Wref  at  t  =  66  s,  the  control  program 
shuts  the  electrolyzer  off.  W(t)  then  starts  to  rise  (from  t  = 
66  s  t o  t  =  135  s).  Figs.  4  and  5  show  that  there  is  a  surplus 


Fig.  6.  Comparison  of  the  electrolyzer  current  simulated  and  measured. 
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Fig.  8.  Battery  energy  W(t). 
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power  on  the  dc  bus  during  this  period,  and  this  surplus  is 
transferred  into  the  battery.  From  t-  136  s  to  t  =  200  s,  W(t) 
decreases  again  due  to  the  high  load  current  (Fig.  5)  and  the 
insufficient  energy  supply  from  the  wind  turbine.  At  t-  295  s, 
W{t)  reaches  the  threshold  value  W\ow  at  which  the  control 
program  turns  the  FC  on  (Fig.  7),  to  supply  the  required  en¬ 
ergy  from  the  stored  hydrogen.  Beyond  t  =  350  s,  the  WT 
current  is  essentially  zero,  while  the  load  current  remains  at 
about  25  A  and  the  FC  current  at  only  17  A,  and  hence  W(t) 
decreases  again  due  to  the  net  8  A  drawn  from  the  battery. 
The  difference  between  the  simulated  and  the  measured  W(t) 
seems  to  be  less  than  2%.  This  is  a  good  indication  of  the 
reliability  of  the  models. 

5.  Residential  energy  consumption  scenario 

A  realistic  scenario  of  the  energy  consumption  of  a  resi¬ 
dence  is  now  developed.  The  values  used  for  the  power  pro¬ 
ductions  of  the  PV  array  and  of  the  wind  turbine  are  typical 
for  our  region.  The  net  dc  bus  power,  that  is  the  algebraic 
sum  of  powers  from  the  PV  array,  the  wind  turbine,  and  the 
load  (this  one  negative)  but  without  the  contribution  of  the 
FC  and  of  the  electrolyzer  is  shown  in  Fig.  9.  It  stays,  on 
the  average,  negative  until  t  =  3800  s,  and  slightly  positive 
afterwards.  Positive  power  values  mean  that  excess  power  is 
available  at  the  dc  bus,  while  negative  power  values  means 
that  there  is  a  power  deficit  at  the  bus.  The  simulation  of  this 
scenario  covers  2  h  (7200  s). 

The  electrolyzer  and  FC  model  parameters  are  again  given 
in  Table  2,  while  their  on/off  control  criteria  are  similar  to  the 
previous  simulation.  The  target  value  Wre f  is  0.70  Wo.  The 
electrolyzer  is  started  when  W(t)  >  0.75  Wmax  and  stopped 
when  W(t)  <  0.70Wmax.  The  FC  is  started  when  W(t)  < 
0.65  Wmax  and  stopped  when  W(t)  >  0.70Wmax.  The  maxi¬ 
mum  capacity  of  the  battery  ( Wmax)  is  set  at  1 1  kWh.  The 
initial  battery  charge  (Wo)  is  set  at  0.75  Wmax. 


Time  (s) 

Fig.  9.  Net  power  on  the  dc  bus  without  the  contributions  of  FC  and  elec¬ 
trolyzer. 


Fig.  10.  Battery  energy  W(t):  Whi  =  0.75  Wmax,  Wref  =  0.70 Wmax,  Wiow  = 
0.65  Wmax. 


The  evolution  of  W(t)  is  displayed  in  Fig.  10.  When  the 
system  starts  (t = 0  s),  W{t)  is  at  0.75  Wmax,  and  the  electrolyzer 
is  running  (Fig.  11).  The  battery  energy  W(t)  decreases  (see 
Fig.  10)  because,  again,  the  net  power  (Fig.  9)  is  negative 
(the  consumption  is  more  than  the  energy  production)  and 
because  the  electrolyzer  is  also  running.  At  t  -  990  s,  W(t) 
reaches  Wref  =  0.7  Wmax  (Fig.  10),  and  the  control  algorithm 
stops  the  electrolyzer  (Fig.  11).  From  t  =  991  s  to  t  =  2995  s, 
neither  the  FC  nor  the  electrolyzer  are  operating  (Fig.  11), 
and  W(t)  keeps  decreasing,  since  the  net  power  on  the  dc  bus 
is  negative  during  this  period  (Fig.  9).  W(t)  keeps  decreasing 
until,  at  t-  2995  s,  it  reaches  the  value  of  W\ow  =  0.65  Wmax, 
at  which  point  the  control  logic  turns  the  FC  on,  and  energy 
flows  back  into  the  battery  from  the  converted  hydrogen.  The 
FC  is  operating  from  then  until  6900  s  (Figs.  10  and  11),  at 
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Fig.  11.  Power  of  the  electrolyzer  and  the  FC. 
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Fig.  12.  Energy  of  the  stored  hydrogen 

which  point  W(t )  reaches  WYQ f  again  and  the  FC  is  turned 
off. 

Fig.  12  shows  the  equivalent  evolution  of  the  energy  of 
the  stored  hydrogen.  The  initial  stored  energy  (Wh?  0)  is  set 
to  15  kWh.  It  is  obvious  that  when  the  electrolyzer  is  in  oper¬ 
ation  producing  hydrogen  (0  <  t  <  991  s),  Wu2(t)  increases, 
and  when  the  FC  works  consuming  hydrogen  ( t  >  2990  s), 
Wh2(0  decreases. 

These  results  show  that  the  models  described  in  this  paper 
can  be  used  to  predict  the  performance  of  a  renewable  energy 
system  with  hydrogen  storage  for  any  reasonable  scenario. 
Moreover  it  is  possible  to  know  the  conditions  necessary  to 
make  the  system  autonomous  by  checking  the  hydrogen  suf¬ 
ficiency  in  the  system  (Fig.  12).  The  stored  amount  of  hydro¬ 
gen  energy  gives  us  the  time  during  which  the  RESHS  could 
work  as  a  purely  stand-alone  source  in  the  extreme  case  of 
zero  power  input  from  the  WT  or  the  PV  array. 

6.  Conclusions 

We  developed  a  model  to  simulate  a  power  generator  fed 
by  renewable  energies  (wind  and  light),  with  batteries  and 
gaseous  hydrogen  as  energy  reservoirs,  and  electrolyzers  and 
fuel  cells  as  converters  of  energy  between  electrical  and  hy¬ 
drogen.  Special  attention  has  been  given  to  the  modeling  of 
each  sub-unit  of  the  system.  The  electrolyzer  and  the  fuel  cell 
are  modeled  by  using  their  steady- state  polarization  curves. 
Due  to  the  important  role  of  power  interfaces  (boost  and  buck 
converters)  in  the  RESHS,  they  are  modeled  for  both  their 
transient  and  steady  state  behaviors  (ARM A  model).  A  sim¬ 


plified  model  is  used  for  the  battery,  with  parameters  easy  to 
estimate,  and  tested  with  all  the  other  system  components. 
The  validation  of  all  component  models  is  based  on  a  10- 
stage  algorithm  that  develops  their  dynamic  evolution  during 
the  course  of  the  simulation.  Comparison  with  experimen¬ 
tal  data  confirms  that  these  models  give  realistic  and  reliable 
results.  The  simulation  gives  an  average  deviation  estimated 
at  less  than  5%  compared  to  an  actual  test-bench  generator, 
even  in  the  presence  of  strong  fluctuations  of  the  primary  en¬ 
ergy  input  sources.  A  second  application  of  the  model  was 
to  simulate  the  power  scenario  for  a  residential  application, 
with  power  sources  typical  for  our  region.  The  results  confirm 
again  that  the  model  can  describe  the  behavior  of  a  RESHS 
in  a  realistic  way. 
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